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The scattering phase shift encodes a good amount of physical information which 
can be used to study resonances from scattering data. Among others, it can be 
used to calculate the continuum density of states and the collision time in a res- 
onant process. Whereas the first information can be employed to examine the 
' evolution of unstable states directly from scattering data, the second one serves as 

l' a tool to detect resonances and their properties. We demonstrate both methods 

' concentrating in the latter case on 'exotic' resonances in tttt and irK scattering. 

Ph. 
& 

1. The Phase Shift and its Ambiguities 

• i— I , 

. The phase shift Si is a convenient variable to parameterize the scattering 

Vh^ I amplitude as known from many textbooks on quantum mechanics. In prin- 

ciple, the phase shift can be extracted from the differential cross section 
da/dfl or other observables. To understand a physical problem from a 
scattering experiment it is often sufficient to know the differential cross 
section. A good example is a bump in the cross section which signals often 
a resonance. However, the knowledge of the phase shift (or amplitude) can 
give us additional information on the dynamics of the scattering process, 
the properties of occurring resonances etc. We shall discuss two such appli- 
cations of the phase shift in the subsequent sections. Both have to do with 
resonances as intermediate states in the scattering process from which Si is 
determined, ft is therefore not unimportant to remind the reader of some 
facts regarding the extraction of the phase shift from experiment. Only if 
the extraction is unambiguous and clean can we rely on the information 
stored in the extracted data. The usage of the differential cross section is 
not enough for this purpose because of several ambiguities encountered in 
the process. Here we will just mention two of them: the so-called Yang 1 
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and Minami 2 ambiguities. The first one allows for two different solutions 
in case that only S and P waves are present (important) in the reaction 
under consideration. The relation between the two solutions, primed and 
unprimcd, reads: S(S 31 ) = 6'(S 31 ) and 8(P 33 ) - 5(P 31 ) = S'(P 31 ) - S'(P 33 ). 
The second famous ambiguity is called Minami after its discoverer. It is 
based as all such ambiguities on the invariance of the differential cross sec- 
tion under a certain transformation. In case of the Minami ambiguity this 
is Sj + i/2 <-> Sj -i/2 where j is the total angular momentum. Hence more 
specifically one obtains 6 (Si) <-» 5 (Pi), S(P 3 ) *-* S(D 3 ) and 6(D$) <-» 5(F§). 
There are more ambiguities which we will not discuss here. But the point 
we wish to make is that because of these ambiguities a differential cross sec- 
tion is not sufficient to extract the phase shift and more data such as the 
polarization observables are necessary 3 . In addition, the extracted phase 
shifts must pass a number of tests based on basic properties of the ^-matrix 
like unitarity, analyticity and crossing symmetry. For instance, in tttt scat- 
tering these are gathered in the so-called Roy equations and the mr phase 
shift has to satisfy the latter in order to be accepted as an unambiguous 
physical quantity 4 . 

This way the experimental phase shifts are indeed trustworthy physical 
values which we can use to get some insight into the scattering process. In 
the next section we will give some examples on how these values might be 
used and processed. 

2. Physics from the Phase Shift 

Having convinced ourselves that the extraction of the phase shift is reli- 
able, we proceed to discuss two examples of how to use and interpret the 
derivative of the phase shift to obtain physically relevant information 
on the resonances occurring in the process. 

2.1. Continuum Level Density 

While calculating correction factors B and C to the equation of states of an 
ideal gas, namely, pV = RT[1 + B /V +C/V 2 }, Beth and Uhlenbeck 5 found 
that the derivative of the phase shift is proportional to the difference of the 
density of states (of the outgoing particles) with and without interaction. 
In case of the I th partial wave we have then, 

m(E) - n\ ' (E) = — (1) 
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This is an interesting result. To appreciate it, let us briefly recall the Fock- 
Krylov method to study the time evolution of unstable states. It is based 
on the fact that unstable states cannot be eigenstates to the Hamiltonian 
and as a result we can expand the resonance states in terms of the energy 
eigenstates; i.e. 

|*) = J dEa(E)\E) . (2) 

With a little bit of algebra we can recast the survival amplitude in a Fourier 
transformation of the so-called spectral function p^ = |a(^4)| 2 which is the 
probability density to find the state \E) between E and E + dE. This way 
we obtain the Fock-Krylov celebrated result 6 

r-OO 

A 9 (t)= / dEp^(E)e- iEt (3) 



used in many investigations on quantum time evolution 7 . Now the proba- 
bility density and the continuum density of states are related by a constant. 
As long as there are no interfering resonances ^ is positive and the above 
identification works without doubt. The best example is the Lorentzian 
(Breit-Wigner) spectral density which is very often used in investigations 
of the decay problem 8 . Indeed, T = & ves Tise t° a phase shift 

S such that ^§ = (E R -Ey+r 2 /4 ' ^ n case 01 severa l resonances the interfer- 
ence pattern can cause the derivative of the phase shift to become negative. 
Therefore, it is safe to take the first resonance and neglect the subsequent 
contribution of the higher lying resonances. We can do that if we want to 
study only the large time behaviour of the time evolution because in this 
case due to the time-energy duality we need to know only the threshold 
behaviour of the phase shift. We can parameterize p^ without loss of gen- 
erality in the form p^ oc (E — -Eth.) 7 to account for the threshold. Hence 
we have 

p* i (E)=G(E)(E-E th .yVK (4) 

The form factor has the property G(E) — > sufficiently fast as E — ► oo and 
a complex pole z , i-e. 1/G(zq) = such that ^sm(z ) < and 5fte(zo) > 0. 
This pole represents the resonance at ER—iT/2 and leads to the exponential 
decay. This generalizes the simple Breit-Wigner form which has the pole but 
no threshold behaviour. These general properties of the spectral function 
allow us to compute the survival probability by going to the complex plane. 
We choose the closed path Cr = Cq + Csr + C^ 4 , starting from zero (after 
change of variables y = E — E t h) along the real axis (Csr) attaching to it 
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a quarter of a circle with radius R (C 1 / 4 ) in the clockwise direction and 
completing the path by going upward the imaginary axis up to zero (Cq). 
In the integral we let R go to infinity noting that along Cjjj{^ , the integral 
is zero. We subtract the contribution along the imaginary axis. This gives 

A 9l (t) = A* l {t) + Al l {t) (5) 

with 

A%(t) = e - iEth t lim / dze- lzt z^G{z + E th ) = C ie - lEnt e- T / 2t (6) 
R ^°° Jc K 

by Cauchy's theorem and 

Al t (t) = C 2 e- iEth t / dxe- xt x^G{-ix + E th .) 
Jo 

cC 2 e- iE ^ t G(E th .)T( 7 + l)-^ (7) 

for the integral along Cq, where the approximation is valid for large times t. 
In the above, r(x) is the Euler's gamma function and C\ is a constant easily 
calculable in terms of the parameters Er, V and E t h. and C 2 is (— i) 7+1 . 
Equations (6) and (7) show the general features which we would expect: an 
exponential decay law followed by inverse power law corrections. The latter 
is independent of the details of G(E) displaying also the dual nature of time 
and energy. Our choice of the elastic resonant reaction from which we want 
to extract the information on the long tail of the decay is dictated by very 
good data at threshold. We opted for a + a — > 8 £>e(2 + ) — > a + a (see 
Figure 1) 9 . The analysis of this data following the mathematical method 
outlined above or alternatively by numerical integration reveals that the 
survival probability behaves as Py(t) = |^4*(t)| 2 ~ t~ 6 - 36 for large times 

10 

Interestingly one can also use the interpretation of dd / dE as continuum 
level density to calculate the density of resonances per unit volume and unit 
invariant mass in a thermal environment e.g. in heavy ion collisions 11 . This 
is m K / iSjh JM cxp(g p 1 /T)±i and is Proportional to the probability 
density to form this resonances. 

2.2. Time delay 

In Figure 1 we can see that all resonances of s Be with I = 2 are nicely 
mapped through the peaks of the derivative of the phase shift and the 
positions of these peaks correspond to masses of the resonant states. We 
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Figure 1. D-wave phase shift (upper half) [9] and its derivative (lower half) in a-a 
elastic scattering as a function of E ex = E c . m . — -E8Be(groundstatc) • The figure on the 
right displays the region of the first 2 + level of 8 Be in detail. The inset displays the 
accuracy of our fit near the threshold energy region which is crucial for the large time 
behaviour of the decay law. 



would expect this due to the interpretation of dS/dE as a continuum level 
density. As explained above, such a level density should have a complex 



given at the International Workshop PENTAQUARK04, July 20-23 at Spring-8, Japan 



6 

pole responsible for the exponential decay at the resonance position which, 
in the experimental data is reflected through a bump in the cross section. If 
we want to map all resonances by this method, then it is more appealing to 
re-interpret dS/dE as a collision time or, in case this is positive, as a time 
delay in a scattering experiment. Such an interpretation was pioneered by 
Wigncr, Eiscnbud and Bohm 12 and is a topic of standard textbooks by 
now. For a wave-packet A(E',E) centered around E the exact expression 
is 13 

f°° dt) 

At(E) = 8n 2 J o dE'\A(E>,E)f2— (8) 

which for a sufficiently narrow wave-packet A(E' , E) gives 

, > dS , „ 

At(E) = 2— . (9) 

With this interpretation we can reinforce the expectation that the collision 
time d5/dE peaks in the vicinity of a resonance (at the resonant energy to 
be exact). Certainly, a collision is delayed if an intermediate state becomes 
on-shell (this happens usually in the s-channel, but some curious examples 
of t-channel singularities also exit 14 ). We emphasize that the collision time 
(9) is strictly the difference between time spent with and without inetraction 
and not simply the time that a projectile spends in the scattering region of 
radius a. The latter in the presence of interaction reads 13 

POO 

T(a) = At(E) + 8tt 2 / dE'\A{E', E)\ 2 
Jo 

and is positive definite in contrast to (9) which we expect to be positive at 
resonance energies. Equation (10) is only applicable if we can define unam- 
biguously a radius a which has a semi-classical character and/or assumes 
that the scattering particles have some extension which, as the example of 
e+e~ — > Z a — > e + e~ shows, is not always a good assumption in resonant 
scattering. Furthermore, resonance production is per se due to interaction 
and therefore (9) the right concept to use for our purposes. 

We can now say that the survival amplitude (3) of a resonance is a 
Fourier transform of a collision time in momentum space if this resonance 
is an intermediate state in the process. 

3. Time delay and resonance physics 

Having identified the derivative of the phase shift as continuum level density 
and as time delay in resonant scattering, we can proceed to apply this con- 
cept to realistic examples (one of them is already displayed in Figure 1). It 



2a - 



sin 2(5 + k'a) 
k' 



(10) 
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is, however, instructive to dwell first on some theoretical connections, mis- 
conceptions and expectations. We note that we consider the usage of time 
delay in resonance physics as a supplementary tool to the other established 
methods. 

In literature one often encounters the statement of the correspondence 
'phase shift motion' <-> resonance. Time delay is nothing else but the exact 
mathematical formulation of this correspondence. However, this correspon- 
dence often carries a misunderstanding as it is attached also to a 7r-jump 
of the phase shift. We stress that this 7r-jump is not a necessary condition 
for a resonance. In the spirit of time delay the condition is a peak around 
the resonance energy. Indeed, there are examples of prominent established 
resonances without the strong 7r-jump like n + a — > 5 _ffe(Piy 2 ) n + a 
which is purely elastic with a jump from 0° to 40° 'only' 15 . 

A simple Breit-Wigner parameterization of the amplitude, i.e. T = 



r,(/w f /■') 



which gives 



e r -e-iT/2 corres P on ds to S = \ tan 
(z§)b-b r = F Br^i/2 ■ ^ ms wom d mean that time delay is negative if 
Br < 1/2. An improvement can be reached by including a non- resonant 
background parametrized here by the diagonal phase Q 16 and energy de- 
pendent width. One then gets {%) e=Er = r( ^ fl) 3^-1/2 + e=Er 
which, in principle, can save the time delay from becoming negative near a 
resonance. However, we would not expect that when the resonant contri- 
bution is large. 

Let us now confront this with experiment. In Figure 2 we have plotted 
the phase shift for the Sn resonances, the inelasticity parameter (note that 
in case there are several channels, the S'-matrix is written as ryexp(2z(5)) 
and the time delay. First of all we find sharp peaks at 1.5 GeV and 1.65 
GeV corresponding to the well known resonances (Particle Data Group 
estimate of the pole value of the first Sn resonance is 1.505). Secondly, 
we get these peaks in spite of the small branching ratio of 5n(1535) which 
is Br(irN) = 35 - 55% and Br(r)N) = 30 - 55%. It is also clear that 
the time delay becomes negative when the inelasticity parameter is largest. 
This can be understood as the loss of flux from the elastic channel due 
to the interpretation of At as density of states 17 . In Figure 3 we have 
done a similar exercise for the Pn case 18 . This is interesting from several 
points of view. Again we find two established resonances, but the focus 
is here on the three star Pn(1710). We find this resonance by the time 
delay method at the right position even if the 7rN branching ratio is as 
small as 10 - 20%. We find it by using the FA02 19 amplitudes even if 
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Figure 2. Single energy values of (a) inelasticities (triangles), (b) phase shift (circles) 
and best fit curves to phase shift (solid lines) in the Sn partial wave of ttN elastic 
scattering and (c) distribution of time delay as a function of energy available in the ttN 
centre of mass system. 



TtN elastic scattering 



FA02 Poles Branching 11 
(GeV) (TtN mode) 




Figure 3. Time delay in the P\\ partial wave of ttN elastic scattering evaluated using 
the FA02 T-matrix solutions. 



the group which has performed the FA02 partial wave analysis cannot find 
the pole corresponding to Pii(1710). Whether or not this resonance exists 
is important for the theoretical prediction of the mass of the Pentaquark 
0(1540) 20 . Indeed, this prediction relies on the existence of the Pn(1710) 
21 . Through the time delay method we find this resonance and also the 
Pentaquark 22 at the right positions. 

In passing we note that even resonances like Pi3(1585), Gi7(2190) and 
if i 9 (2220) with ttN branching ratios of 10 — 20% leave clear fingerprints in 
the time delay plots 18 . 
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4. Resonances in tvtv and tvK scattering 

The previous sections showed that the time delay method is reliable in 
nuclear and baryon resonance physics. We now turn our attention to the 
mesonic case 23 . To show how reliable the method indeed is and how sen- 
sitive it is to small phase shift motion, we first apply the method to the 
case of the p-mesons. This is depicted in Figure 4. Evidently, we find the 



200 




Figure 4. Time delay plot of the p resonances evaluated from a fit to the p-wave phase 
shifts in tttt elastic scattering. 



'not-to-be-missed' p(770), its first excitation p(1450) and its second exci- 
tation p(1700) which are all indisputable resonances. The peak at 1166 
McV corresponds to a small phase motion and one could be tempted to 
disregard it as a fluctuation. However, several other cases, among oth- 
ers the three star resonance £>i3(1700) and the two star Fi 5 (2000), show 
that small phase shift 'motion' can signify a resonance. This seems to be 
the case also here. Particle Data Group lists also several mesons between 
1100 — 2200 MeV which by itself is not a remarkable fact. But at the recent 
Hadron 2001 conference in Protvino some authors have pointed out a grow- 
ing evidence for a p-like resonance at 1200 MeV which we think appears in 
our time delay plot . Such a state has been also predicted in a coupled 
channel quark/meson model 25 . Our result in the p- meson sector is then an 
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independent confirmation by the time delay method. 

In the last few years the scalar sector attracted lots of attention. One 
of the reasons is the 're-discovery' of the famous cr-meson and its 're- 
appearance' in the Particle Data Book. The difficulty with this meson is 
reflected in the wide range of its possible mass, 400 — 1200 MeV. The time 
delay analysis for this sector is summarized in Figure 5. Of course, /o(980) 
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Figure 5. Time delay plot of the scalar meson resonances evaluated using two different 
sets of the s-wave phase shifts [26,27] in 7™ elastic scattering. 

is a dominant contribution here. We identify the peak around 1.23(1.34) 
GeV with / (1370) for which Particle Data Group quotes the range of pos- 
sible pole mass between 1.2 and 1.5 GeV. Similarly the peak at 1.42(1.47) 
GeV is attributed to / (1500) (the PDG value is 1.4-1.6 GeV). The anal- 
ysis of both phase shifts reveals a resonance at 700 MeV. If, in addition, 
we take the information of the Kaminski phase shift we see also a peak at 
350 MeV. Can we take this as an evidence for two resonances? Let us first 
note that in the region of 400 — 1200 where the cr-meson is found, one can 
identify two accumulation points. One at 350 — 600 MeV and the other one 
at 700 — 850 MeV. The low lying case is supported also by unitarized meson 
model 28 , unitarized chiral perturbation theory 29 and by unitarized quark 
model 30 , by coupled channel analysis 31 and by the so-called ABC effect 
32 which is with us since 1961 and the recent decay J/ty — > <toj — > 7r7ra; 33 . 
The 700 MeV case finds its confirmation in Nambu-Jona-Lasinio models 34 , 
Weinberg's mended symmetry 35 and Bethe-Salpeter calculation 36 . Hence, 
these two accumulation points are not artificial constructs. They can be 
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seen from the experimental results the theoretical expectation 37 so from the 
time delay method. Figure 6 serves as a cross check if the resonance found 
by the time delay method in the elastic channel finds its confirmation also 
in other channels. As can be seen this is indeed the case. 

In the strange scalar sector the controversy regarding the lightest scalar 
(called k meson) is even bigger. Figure 7 is an analysis of this sector using 




1.1 1.2 1.3 1.4 1.5 1 1.1 1.2 1.3 1.4 

E (OeV) E (GeV) 



Figure 6. Energy dependence of the s-wave phase shifts and time delay in the (a) 
■kit — > KK and (b) K K —* KK reactions. The phase shifts (solid lines in (b)) have 

been determined by subtracting the fit to the CERN-b rnr phase shift data from <j> n K i n 
the above figure (a). Dashed lines show the same quantities in the model calculation of 
[40]. 

two different phase shifts. This analysis also reveals the existence of two 
low lying resonances: one at 0.8 GeV and the other around 1 GeV which 
we identify with the putative K-meson. 

We have applied the time delay method to many 'standard' cases, the 
established baryon resonances, the p mesons and the K*- as well as K%- 
mesons (discussed in 23 ) and found a good agreement with data. Some less 
established resonances found by different methods get confirmed through 
the time delay method. By using the KN phase shift we found the recently 
discovered Pcntaquark with a mass very close to the observed and predicted 
value 22 . We found the spin-orbit partners of this Pentaquark very close to 
the theoretical expectations 38 , 39 . Last but not least, our nuclear physics 
case discussed here in section two, shows also the virtues of the time delay 
method not only in finding nuclear levels, but also in studying the quantum 
evolution of unstable systems for large times. 
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Figure 7. Time delay plots evaluated from the phase shifts in s-wave, isospin 1/2 Kn 
scattering, displaying the strange scalar mesons. 
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